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CHIP-FIRING ON GENERAL INVERTIBLE MATRICES 


JOHNNY GUZMAn AND CAROLINE KLIVANS 


Abstract. We propose a generalization of the graphical chip-firing model 
allowing for the redistribution dynamics to be governed by any invertible in¬ 
teger matrix while maintaining the long term critical, superstable, and energy 
minimizing behavior of the classical model. 


1. Introduction 

The classical chip-firing model describes the dynamics of a Laplacian action on a 
graph and can be seen as a discrete diffusion process, see e.g. [smui]. For a graph 
with n vertices, configurations are integer vectors c S Z” where Ci is interpreted 
as the number of chips at vertex i. At each time step, chips are redistributed via 
local rules. For a configuration c, if a vertex i has at least as many chips as it 
has neighbors, i fires sending one chip to each of its neighbors, forming the new 
configuration c'. This action is easily identified as subtracting the appropriate row 
of the graph Laplacian: c' = c — L^ e^, where L is the graph Laplacian and is 
the ith indicator vector. 

This model is well studied and has a strong theory of long term behavior in terms 
of self-organized criticality [Si, superstability [mili] and energy minimization [U 
I12j . In particular, every initial configuration eventually stabilizes to a unique stable 
configuration - the order of firings does not matter. 

Already in m, Gabrielov broadened the setting for chip-firing realizing that one 
could consider a finite set of states and the dynamics governed by any avalanche 
finite matrix. The system did not need to come from an underlying graph; the 
firing rule is governed by the matrix. Such systems similarly converge to a set of 
stable configurations regardless of the order of firings. 

In more recent work, the present authors |12j recognized avalanche-finite ma¬ 
trices as precisely M-matrices, a well known class of matrices from other contexts 
and continued to expand the theory of chip-firing in this setting. In particular, 
the theory of criticality, superstability and energy minimization extends to all M- 
matrices. The name avalanche finite, however, is quite appropriate. If one attempts 
to chip-fire using a matrix outside of this class, Gabrielov showed that not all initial 
configurations will eventually stabilize. Some state will always be able to fire. 

In this paper we propose a significant broadening of the chip-firing model allow¬ 
ing for the redistribution dynamics to be governed by any invertible integer matrix 
while maintaining the long term critical, superstable, and energy minimizing behav¬ 
ior of the classical model. To overcome the seeming contradiction with Gabrielov’s 
work, we change the scope of valid configurations over which chip-firing occurs. In 
classical chip-firing, valid configurations are precisely non-negative integer vectors, 

Date: August 19, 2015. 

Key words and phrases, chip-firing, energy minimization, combinatorial Laplacians. 


1 



2 


J. GUZMAN AND CAROLINE KLIVANS 


i.e. integer points in the closure of the positive orthant. In our more general setting, 
the positive orthant is replaced by a different cone. As a consequence, valid (and 
hence also critical and superstable) configurations may contain negative entries. 

The general set-up is as follows. For a fixed invertible matrix L, one must 
choose an M-matrix M. The chip-firing model will be dictated by the pair (L,M). 
The matrix M can be any M-matrix, it need not depend on L. Hence for a fixed 
distribution matrix L, one can in fact define many chip-firing models. The valid 
configurations are a collection of integer points, denoted S'^, determined by the 
matrix LM“^. Given a configuration c £ S'^, a site i is allowed to fire if the 
resulting configuration c' = c — Le^ remains in S~^. The well-behaved dynamics 
of this model are established by relating chip-firing via L over S~^ to chip-firing 
via M over a related set of (not necessarily integer) configurations. As in the 
graphical case, critical, superstable and energy minimizing configurations of all 
form systems of representatives of the critical group of L, K(L) := coker(L). 

As a special case, one easily recovers classical chip-firing on a graph. Let L be 
(the transpose of) a reduced graph Laplacian. Then L is itself an M-matrix. Thus 
for such an L (or for any initial M-matrix) one can choose M = L giving the pairing 
(L,L). As we will see in Section 


6.1 


S~^ then reduces to Z”p and we recover the 
well known chip-firing model for graphs and digraphs. 

As another special case, one can always set M = I, since the identity matrix is 
an M-matrix. The pairing (L,I) is particularly interesting. The superstable config¬ 
urations for (L, I) are precisely the integer points of the fundamental parallelepiped 
of the lattice generated by L, see Section [6?^ It is well-known that these integer 
points form a system of representatives for coker(L). The theory here now gives a 
dynamical system which achieves these configurations. 

In Section we explain the general set-up for chip-firing over an invertible 
integer matrix L with respect to a choice of M-matrix M. In Sections[^|^ andj^we 
show the existence and uniqueness of critical, energy minimizing, and superstable 
configurations respectively. Finally, in Section we consider how the long term 
dynamics vary under different choices of M including M = L and M = I. We also 
consider the special case of chip-firing in higher dimensions where L is a reduced 
combinatorial Laplacian of a simplicial complex as in [10] . 


2. General Chip-Firing 

Let L be an n X n invertible matrix with integer coefficients. Consider n as the 
number of states of a system. A configuration of the system is any integer vector 
c S Z". The value is the amount of commodity (chips, dollars, flow, etc.) at site 
i. L dictates the redistribution of the commodity among the sites of the system. 

Now let M be an n X R M-matrix with real entries and define N := LM“^. The 
matrix N will dictate which configurations are valid configurations of the system. 
In particular, given N define the set: 


S+ = {NxiNa; e Z”,a: € Mmol¬ 
s'’'" consists of all integer vectors formed by multiplying N by a non-negative real 
vector. Note that without the non-negativity condition on x, this set would simply 
be all of Z” since N is invertible. The set of integer points in S~^ are the valid 
configurations for the pairing (L,M), under the dynamics of L. In classical chip¬ 
firing on a graph, S'^ is exactly the non-negative integer points. 
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For a valid configuration c G , a site i is ready to fire if the configuration 
c' = c — Lei, resulting from subtracting the ith row of L, is also valid, i.e. if c' 
remains in 5'“'". A conhguration is called stable if no site is ready to fire. Similarly, a 
set of sites can multihre (with multiplicity) if the resulting conhguration c" = c—Lz 
is also valid, i.e. if c" remains in S'"*", where Zi is the number of times site i hres. A 
conhguration in which no set of sites can multihre is called superstable. We will in¬ 
vestigate the properties of superstable conhgurations in the following sections along 
with notions of critical and energy minimizing conhgurations. In particular, we will 
consider these conhgurations with respect to the following equivalence relation on 
conhgurations. 

Definition 2.1. For two conhgurations f,gG Z", dehne / if / = 5 — Lz for 
some z e Z". 


Namely, the equivalence classes form the elements of the cokernel of L, coker(L) 
= Z^/imL, the critical group of L. Each of the collections: superstable, critical, 
and energy minimizers will form a system of representatives under this equivalence 
relation. Namely, there exist unique critical, energy minimizing and superstable 


conhgurations per equivalence class in S'^, see Theorems 3.5 4.3[ and 5.5 


These results will follow naturally from a second chip-hring system running in 
parallel with distribution matrix M. Importantly, M is an M-matrix, which are 
widely used in many contexts such as economics and scientihc computing, see e.g. [6l 
HIS] and m and references therein. There are many known characterizations of 
M-matrices. We mention the three most important in our context. 


Definition 2.2. An n x n non-singular matrix M such that Mij < 0 for all i ^ j 
and Mii > 0 for all i, is an M-matrix if any of the following equivalent conditions 
hold: 

(a) M is avalanche hnite. 

(b) All entries of are non-negative. 

(c) There exists x G K” such that x > 0 and Mt has all positive entries. 

The equivalence of condition (a) is due to Gabrielov [TT] and refers to a system 
with non-negative valid conhgurations. We refer to Plemmons m for conditions 
(b) and (c) and many others. 

The valid conhgurations corresponding to a choice of M are the vectors that 
generate points in S'+. Specihcally, dehne 

R+ = {xGE.'^\NxGZ'^,x>0}. 

The elements of i?'*" are non-negative real vectors such that Na; is integer, equiv¬ 
alently Na; is in 5'+. The elements of are the valid conhgurations for the pair 
(L, M) under the dynamics of M. 

Definition 2.3. For two conhgurations x,y G K", dehne x ~m y x = y — Mz for 
some z € Z". 


Proposition 2.4. Let f,g G S'^ and x,y G such that f = Na; and g = Ny. 
Then 


f g ijf x-^M y- 
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Proof. Suppose that / ~l 9 then we have the following relations: 

f = g-Lz 
Nx = Ny — Lz 
N"^ Nx = N"^ Ny - N"^ Lz 
X = y — ML“^ Lz 
X = y — Mz. 

□ 


Hence we see that (multi-)firing on configurations in S'^ via L is equivalent to 
(mulit-)firing on conhgurations in via M. The chip-hring dynamics in the second 
case are not the same as chip-hring on M-matrices as in [12] , where the set of valid 
conhgurations is always Z>q. However, we will see that much of the nice behavior 
carries over when considering a general . We end this section with an example. 


Example 2.5. Let 


2 

-1 


/ 3 

-1 


/ 1 

.25 

1 

2 

—1 and M = 

-1 

3 

-1 

. Then N = -.25 

.5 

1 

-1 

2/ 

1-1 

-1 

3/ 

1 .75 

.25 


.75> 


The valid conhgurations, S~^, consist of integer points of the form Nx such that 


X > 0. Vectors in 5'+ include 


On the other hand, 


and 


N- 



-1.75 

.25 ] ^0 and N"^ 
2.25y 



3. Criticality 

We hrst investigate the long term behavior of successively bring one state at a 
time. Assume throughout that the total number of sites in the system is n and let 
Ci be the ith standard basis vector of K". A site is stable if no individual site can 
validly hre: 

Definition 3.1. A conhguration / G S'"*" is stable if / — Le^ ^ S'^ for 1 < z < n. 

Definition 3.2. A conhguration x G R~^ is stable if x — Mci ^ R'^ for 1 < z < n. 

Note that if / = Nx with / G S~^ and x G i?+, then / is stable if and only if x 
is stable. 

Definition 3.3. A conhguration / G is reachable if there exists a conhguration 
g G S~^ such that 

(i) g — Lci G for 1 < z < rz and 

(ii) f = g - Y!j=i Lcij such that g - Le^, G S+ for all I < k. 
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Namely, a configuration / is reachable if (i) there exists a sufficiently large con¬ 
figuration g such that (ii) / can be reached from 5 by a sequence of valid individual 
firings. A configuration in i?+ is reachable if there exists a configuration y € 
that satisfies the analogous conditions (i) and (ii) under the M operator. 

Definition 3.4. A configuration in or i?+ is critical if it is both stable and 
reachable. 


Note that if / = Ncc with / G S'"*" and x S i?'*', then / is critical if and only if x 
is critical. 

The original work arising from the sandpile literature focused on critical configu¬ 
rations and the phenomenon of self-organized criticality. Starting at any sufficiently 
large initial configuration and performing valid individual firings, a stable conhg- 
uration is eventually obtained, i.e. there are no infinite sequences of valid firings 
(condition (a) from Definition 2.2). Moreover, from a fixed initial configuration, 
the same stable conhguration is always reached regardless of the choice and or¬ 
dering of the individual firings. Equivalently, critical configurations exist and are 
unique per equivalence class under the equivalence relation defined by the image 
of the toppling matrix. This fundamental behavior was considered by Bak, Tang 
and Wiesenfeld [5] for grid graph networks, by Dhar [9] and Speer m for general 
graphical networks and Gabrielov m for M-matrix systems. The behavior contin¬ 
ues to hold in our framework and follows from checking that the results of mm 
extend to the case of restricted valid configurations. 


Theorem 3.5 (See [TTl Theorems 1.2, 2.9] ). Let (L,M) be a fixed pairing of an 
invertible integer matrix L and M.-matrix M. Then, in (resp. R^), critical 
configurations exist and are unique per equivalence class under the relation 
(resp. --m/ 

Proof Sketch. Suppose that x S Then two conditions hold, Nx S Z” and 
X > 0. Firing a state i yields the new configuration x' = x — Me^. Importantly, the 
new configuration x' is such that Nx' remains integer. To see this, expand Nx' as 
follows: 


Nx' = N(x — Mci) 

= LM"^x - LM"^Me, 
= Nx — Le^. 


Nx is integer by assumption and Lci is integer because L is an integer matrix. Hence, 
starting at a configuration in i?'*', checking to see if the resulting configuration is 
valid only requires checking if it remains non-negative. This is the only condition 
used on the configurations in the referenced Theorems. □ 

Note, this implies that the number of critical configurations is equal to the det(L), 
the size of the coker(L). 


Example 3.6. Returning to the pairing (L,M) from Example 2.5 
configurations (of 5'+) are: 


the critical 
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4. Energy Minimization 


In this section we define an energy form on configurations and consider the rela¬ 
tionship between chip-firing and energy as well as the class of valid configurations 
which minimize energy. These ideas were first developed by Baker and Shokrieh 
in [1] and shown to generalize to M-matrices in [12]. Again, the difference in this 
context is that the valid configurations are the points of and that is no 
longer necessarily equal to Z"q. The technical results of this section will allow us 
to establish the existence and uniqueness of superstable configurations in the next 
section. The presentation follows m- 

Given an M-matrix M and a vector q S R ^, define the following energy, 

E{q) = \\M-W\l (1) 

Below we prove that E(q) has a unique minimizer per equivalence class, see The¬ 
orem |4^ To prove this claim, we relate energy to chip-firing in the following two 
Lemmas. First we make a definition. Given z S Z" define z~^ G Z”q by = Zi if 
Zi > 0 and 0 otherwise. 

Lemma 4.1. If x,y G R'^ and y = x — Mz then the configuration h = x — Mz"*" is 
also in R'^. 


Proof. To show that h is in i?+, we need to show that h > 0 and N/i G TP. The 
non-negativity of h follows verbatim from m Lemma 3.2]. Integrality follows as 
in Theorem 13^ □ 


The next result appears as Lemma 3.3 in |12j . It does not consider the domain 
of the two configurations in question and holds in our more general set-up. 

Lemma 4.2. [T^l Lemma 3.3] Let M be an M-matrix and y = x — Mz, then 

E{y) = E{x) + z*z - 2z* a; = E{x) - z‘z - 2z‘ y. 

We can now prove the uniqueness of energy minimizers within the valid config¬ 
urations of an equivalence class under the relation defined by an M-matrix M. 

Theorem 4.3. Let x be a configuration in R'^, then the minimization problem 

( 2 ) 


min 

y~MX, yGR+ 


has a unique minimizer. 

Proof. Suppose w,v G R~^ and w ~m x and v x both minimize Equatio n 
Then w ~m v and there exists z such that w = u —Mz. Since w is in R'^, Lemma 4.1 
gives that h = v — Mz"*" is also in i?+. Now, Lemma |4^ gives: 

E{h) = E{v) - iz+Y{z+) - 2(z+)‘ M-i h. 

The last quantity, 2(z+)‘ M“^ h, is positive. To see this, note that h G R'^ implies 
h > 0, M~^ is non-negative matrix, and z"*" is an all non-negative vector. Therefore, 

E{h) < E{v) - (z+)‘(z+). 

On the other hand, v is an energy minimizer, therefore z"*" = 0. Now consider the 
relation between w and v, the two supposed distinct minimizers, 

E{w) = E{v) + z*z - 2z‘ M-i V. 

If z < 0, then z*z is positive, 2z‘ M“^ v is negative, and E{v) < E{w) unless z = 0. 
Thus z = 0 and w = v. □ 
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Corollary 4.4. Let f be a configuration in S~^ , then the minimization problem 

min ||L-^g||^ (3) 

g^f,g&S+ 

has a unique minimizer. 

Proof. Suppose h,k G 5'+ and h / and k / both minimize Equationj^ Since 
h,k G S~^ , h = Nx and k = Ni/ for some x and y in . Furthermore, h k 
implies that x V- 

First note that if h and k are both minimizers, then x and y are both minimizers: 
If X is not a minimizer then there exists an m which is a minimizer such that 
m = X — Mz. Multiplying by N gives Nm = Nx — NMz g = h — Lz, where 
g = Nm G But then R(g) < E(h) contradicting the minimality of h. 

On the other hand, ||L“^/i|j 2 = ||L“^fc|j 2 and therefore ||L“^Nx ||2 = ||L“^Ni/|| 2 , 

||M-ix||2 = ||M-iy|ir 

But Equation has a unique minimizer by Theorem 


4.3 


□ 


5. Superstability 


Recall that a stable configuration c is a valid configuration in which no (individ¬ 
ual) site can fire such that the resulting configuration is still valid, i.e. c' = Le^ ^ S'^ 
for all i. Superstability extends this notion to account for subsets of sites firing si¬ 
multaneously and with multiplicity. 

Definition 5.1. A configuration / G is superstable if f — Lz ^ 5+ for every 
z e Z”, z > 0 and z 0. 


Definition 5.2. A configuration x G R'^ is superstable if x — Mz ^ i?+ for every 
z S Z”, z > 0 and z 

Note that if / = Nx with / G S~^ and x G R~'~, then / is superstable if and only 
if X is superstable. 


Example 5.3. Let L be the reduced graph Laplacian of the graph in Figure 
where we have deleted the row and column of the Laplacian corresponding to the 
sink vertex. Let M = L. This choice for M is valid because reduced graph Lapla- 
cians are themselves M-matrices, see Section 6.1 Under the pairing (L, L), the valid 
configurations are all non-negative vectors, = Z>q. The configuration shown 
in Figure is stable - no individual site can fire with the resulting configuration 
remaining valid. The configuration is however not superstable - the two sites con¬ 
taining 2 chips can simultaneously fire with the resulting configuration remaining 
valid. 


Remark 5.4. In the graphical case, if the graph is an Eulerian directed graph, 
then it is sufficient to consider subset firing without multiplicity; i.e., in the defini¬ 
tion above, Z" may be replaced by {0,1}". The more general definition above is 
necessary for any non-Eulerian graph or more generally for any M-matrix, see m 
Section 4.1]. 

As with critical and energy minimizing configurations, superstable configurations 
exist and are unique per equivalence class of the coker(L). We establish this by 
showing that x G R'^ is superstable if and only if it is a minimizer of and 
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Sink 


Figure 1. A stable but not superstable configuration. 


therefore / S is superstable if and only if it is a minimizer of The following 
appears essentially as Theorem 4.6 of m, we reproduce it here for completeness. 

Theorem 5.5. Fix an invertible n x n matrix L and an n x n IsA-matrix M. A 
configuration x S is superstable with respect to M if and only if x is an energy- 
minimizer, i.e., a minimizer of Equation^ 

Proof. Let x be superstable. Let y G be equival ent t o x. Then x = y — Mz for 
some z and h = y — Mz’*' is also in by Lemma 4.1 But x is superstable and 
hence z~^ must be 0. Now we can compare the energies of x and y: 

E{y) = E{x) - z*z - 2z‘ w. 

This implies E{y) < E(x) which is a contradiction unless z = 0. 

Now let X be an energy-minimizer. Suppose x is not superstable. Then there 
exists zq such that x — Mzq G with zq > 0. Let y = x — Mzq G R'^ . Comparing 
energies gives: 

E{y) = E{x) - zqZq - 2zlM~^y. 

This implies E{y) < E(x) which is again a contradiction. □ 


Example 5.6. Returning to the pairing (L,M) from Example 2.5 the superstable 
configurations (of S'+) are: 



In the classical model, if L is the distribution matrix, define to be the vector 
given by D\‘ = Ijh — 1 for all i. Then a configuration c is critical if and only if 
— c is superstable. 

This duality does not hold in the more general case. Moreover, there seems to 
be no simple linear relation. Consider, for example, the critical and superstable 
configurations from Example 3.6 and 5.6 Note that in S'^ there is a component¬ 
wise maximum critical configuration, but subtracting in this configuration does not 
give the collection of superstables. 


Question 5.7. Is there a duality between critical and superstable configurations? 


6. Special Cases 

In this section, we consider three important special cases: (1) L is a reduced 
graph Laplacian (2) M is the identity matrix and (3) L is a reduced combinatorial 
Laplacian of a simplicial complex. 
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6.1. Graph Laplacians (L, L). In classical chip-firing, for example as in [3l|4l[9], 
there is an underlying graph to the system. States correspond to vertices of the 
graph and chips are associated to vertices. The distribution dynamics are given by 
the graph Laplacian D — A, where D is the diagonal matrix of degrees of vertices 
and A is the adjacency matrix of the graph. For a graph with n vertices, the graph 
Laplacian is an n x n singular matrix. To fit our context, one must introduce a 
sink vertex - any vertex can be declared the sink and then is no longer allowed 
to fire. The distribution matrix is a reduced graph Laplacian, the row and column 
corresponding to the sink is deleted. The valid configurations consist of all non¬ 
negative integer vectors. 

This special case is recovered by first noting that reduced graph Laplacians are 
in fact M-matrices mi [II] and then pairing the reduced Laplacian with itself. Let 
L be a reduced graph Laplacian (or more generally any M-matrix). In this case, 
one can choose M = L for the pairing (L,L). Under this pairing, 

N = LL-i=I, S'+=Z^-\ R+=Zl-^ 


and we recover the usual graphical chip-firing dynamics. The stable and superstable 
conditions of not leaving S~^ reduce to the conditions that any (multi-) firing should 
yield a non-negative configuration. Example 5.3 uses the pairing (L,L) on a small 
graph. 


6.2. Fundamental parallelepipeds (L,I). The identity matrix I is an M-matrix. 
Therefore, for any fixed L, we can always take the pairing (L, I). In fact, any positive 
diagonal matrix D is an M-matrix. 

Proposition 6.1. Suppose D is a non-negative diagonal matrix. If M and DM 
are M.-matrices, then the superstable and critical configurations of S~^ are the same 
under the pairing (L,M) and (L,DM). 

Proof. We will show that the collection of valid configurations is the same under 
the two pairings. Let correspond to (L, M) and correspond to (L,DM). 

Suppose / e 5'+, then / = Na; = LM“^a: for some a: > 0. For / to be in S'^, we 
must have / = L(DM)“^?/ for some y > 0. Solving for y, we see that: 

LM"^D"^y = LM"^x 
D~^y = X 
y = Dx. 

Since x > 0 and D is positive, Dx is non-negative and / S 5^. 

Similarly, suppose g G 5"^ then g = LM~^D~^ry for some w > 0. In order for g 
to be in the vector v = must be non-negative, which again holds because 

D is a positive diagonal matrix. □ 

By Proposition |6.l[ the critical and superstable configurations of S'^ are the 
same for (L,I) and (L,D). For the remainder of the section we assume that M = I 
for simplicity. Suppose that M = I, then 

N = L and S~^ = {Lx | Lx S Z", x G R>o}- 

Given an invertible matrix L, define the lattice generated by the columns of L: 

C{L) = {Lx|x G Z"}. 
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The fundamental parallelepiped of C with respect to L is: 

iP(£(L)) = {Lx I a; G K”, 0 < a:, < 1}. 

It is well known that the integer points of the fundamental parallelepiped of a full- 
rank lattice form a system of representatives for Z"/imL, see e.g. [H Chapter 9]. 
The next proposition shows that these points can be seen as stable configurations 
in a chip-firing model. 

Proposition 6.2. Under the pairing (L,I), the critical and superstable configura¬ 
tions of S'^ coincide and are precisely the integer points of the fundamental paral¬ 
lelepiped V{C(L)). 

Proof. Let / G S’"*" be a superstable configuration and let / = Lx. Then x is a 
superstable configuration in i?+. Hence x > 0. Suppose Xi > 1 for some i. Then, 
X — Mci = X — Ci > 0, contradicting the superstability of x. 

Therefore every superstable configuration / is contained in V{C(L)). Further¬ 
more, the number of superstable configurations is the same as the number of integer 
points in 'P(£(L)) since both form a system of representatives for Z"/imL. Hence 
the two collections must be the same. 

The argument above actually shows that any stable point of S'^ must be con¬ 
tained in P(£(L)). As critical configurations are also stable and a system of repre¬ 
sentatives for Z^/imL, they indeed must coincide with the superstables. □ 

Example 6.3. Let L be as in Example |2.5| and let M = 1. Then the superstable 
and critical configurations are: 



6.3. Combinatorial Laplacians. The classical framework for chip-firing is on a 
graphical network with distribution given by the graph Laplacian as discussed in 
Section |6.1| In cni, a model for chip-firing in higher dimensions was introduced. 
The underlying configuration of the system is a simplicial complex and the com¬ 
modities (chips, flow) are associated to ridges (codimension-one faces) of the com¬ 
plex. The distribution dynamics are given by a reduced combinatorial Laplacian 
relating codimension-one faces to maximal faces. The algebraic aspects of chip¬ 
firing, i.e. the form and size of the critical group were worked out in m, including 
the connection to higher dimensional spanning trees of the complex. It remained 
an open problem to give a notion of critical or superstable configurations in this 
context, namely sets of representative for the critical group that reflect long-term 
stability dynamics of the system. 

The paradigm introduced here provides well motivated collections of critical and 
superstable configurations. We illustrate with an example and refer the reader 
to m for a more thorough discussion of chip-firing in higher dimensions. 

Example 6.4. Let A be the two-dimensional boundary complex of the tetrahedron 
with maximal faces: (123,124,134, 234}. Faces of the complex will be oriented, we 
will use the standard orientation induced by the vertices. In our example, this 
directs the edges from smaller vertex to larger vertex. The two-dimensional faces 
are also oriented from smaller to larger vertices. Envision traversing the triangular 
face 124 by starting at vertex 1 moving to vertex 2 then vertex 4 and returning 
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back to vertex 1. Traversing the face in this way, we consider the edges 12 and 24 
oriented in the same direction with respect to 124 and oriented oppositely to the 
edge 14 with respect to 124. The consistency of orientation of two edges contained 
in a common triangle will be reflected by a +1 or —1 in the combinatorial Laplacian. 


2 



Figure 2. The complex A is the boundary of the tetrahedron. 
The three dashed edges represent the sink of the system. 


As shown in m, the sink in this context consists of an entire spanning tree 
of the one-dimensional skeleton of A. Choose the sink to be the one-dimensional 
tree T with edges {12,13,14}. With this choice, there will be three states in the 
system corresponding to the three edges not in the tree: {23, 24, 34}. The reduced 
combinatorial Laplacian is formed by deleting the rows and columns corresponding 
to the edges in the spanning tree. In this case, the reduced combinatorial Laplacian 
L of A with respect to T is the matrix from Example |2.5[ our running example: 



23 

24 

34 

23 / 

' 2 

-1 

1 

24 

-1 

2 

-1 

34 \ 

1 

-1 

2 


Configurations in S'+ are elements of which associate an integer value to 
edges. In this context, it is natural to consider these values as flows along the 
edges (as opposed to numbers of chips). Firing an edge consists of subtracting 
the appropriate row of the combinatorial Laplacian. This can be interpreted as 
diverting flow from an edge across triangles that contain it. 

Computing critical and superstable configurations requires a choice of M. Sup¬ 


pose M = 1. Then we find the configurations of Example 6.3 



Figure |6.4| illustrates the three non-zero critical configurations. The third image 
interprets the negative value in the critical configuration not as a negative flow but 
as flow in the opposite direction of the orientation of the edge. 

One narrative for this system is as a model of traffic flow. The edges represent 
streets and the vertices are intersections. Sink edges represent major boulevards 
that can handle large amounts of traffic. As local streets become congested, traffic 
diverts to neighboring streets. Informally, the system stabilizes when as much 
traffic as possible has been diverted to the major boulevards. The third image of 
Figure 6.4 is a solution in which the directionality of street 24 has been reversed. 


We end with a general question given the new framework for chip-firing. 
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Figure 3. The three non-zero critical and superstable configura¬ 
tions of the tetrahedron with sink equal to all edges containing the 
vertex 1 and under the pairing (L,I). 


Question 6.5. Given L, are there natural choices for M? 


We have seen that if L is an M-matrix, then M = L is a natural choice that yields 
the classical model with non-negative valid configurations. We have also seen that 
we can always set M = I and achieve the integer points of the fundamental paral¬ 
lelepiped. What are natural choices if L is the reduced combinatorial Laplacian of 


a simplicial complex, as in Example 6.4? The matrix M from Example |2.5| is the 
reduced graph Laplacian of the facet-to-facet dual graph of A. But, in general, this 
dual graph will not have the correct size to match L. 
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